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Abstract : 

In this paper, we analyze the effect of geometrical constraint on the 
conformational properties of an infinitely long linear semiflexible polymer 
chain confined in-between two constraints under good solvent condition in 
two dimensions. The constraints are two impenetrable stair shaped surface 
and for two dimensional space, the surface is a one dimensional line. The 
semiflexibility of the chain is accounted by introducing a Boltzmann weight 
of bending energy required to produce each turn in the chain and good 
solvent condition was accounted by using self avoiding walk model of the 
chain. We have calculated exact critical value of step fugacity required 
for polymerization of an infinitely long polymer chain confined in between 
the constraints for different values of separation between the constraints for 
directed version of the model. We have also calculated possible maximum, 
minimum values of the persistent length for such chains and the maximum 
value of bending energy required for each turn in the chain for few values of 
separation between the constraints. 
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1 Introduction 

Biomolecules (DNA &i proteins) live in the crowded, constrained regime. 
Such molecules are soft object and therefore can be easily squeezed into the 
spaces that are much smaller than the natural size of the molecule in the 
bulk. For instance, actin filaments in eukaryotic cell [1] or protein encapsu- 
lated in Ecoli [2] found in nature are the examples of confined biomolecules 
that serves as the basis for understanding numerous phenomenon observed 
in the polymer technology, bio-technology and many other molecular pro- 
cesses occurring in the living cells. The conformational properties of a single 
biomolecule have attracted considerable attention in the recent years due to 



developments in the single molecule based observations and experiments 
[3, 4, 5, 6, 7, 8]. Under confined geometrical condition, the excluded volume 
effect and effect of the geometrical constraint compete with entropy of the 
molecule. Therefore, geometrical constraint can modify the conformational 
properties of the molecules. 

The conformational properties of a linear flexible polymer molecule under 
good solvent condition, confined to flat parallel walls (slit) have been studied 
for past few years; for instance, see, [9, 10, 11, 12, 13, 14, 15] and references 
quoted therein. Whittington and his coworkers [12, 13, 14, 15] used self 
avoiding walk model to study behaviour of a surface interacting flexible 
polymer chain confined between two parallel walls on a square lattice. They 
calculated phase diagram of the polymer chain having attractive interaction 
with the walls by solving the directed walk model exactly [12]. Rensburg 
et al. [15] through numerical calculation using isotropic self avoiding walk 
model showed that phase diagram obtained for a surface interacting linear 
polymer chain confined in between two parallel walls has qualitatively similar 
phase diagram to that obtained by Brak et al. [12] for directed walk model 
of the problem. 

Theoretical understanding of a semiflexible polymer chain confined in a 
cylindrical pore, rectangular shaped micro-channel, tube like narrow chan- 
nel or other kind of confined regime is also discussed in the literature to 
analyze effect of the confinement on the conformational behaviour of the 
chain; see, for instance [16, 17, 18, 19, 20, 21, 22, 23] and references therein. 
However, in present investigation, we have considered a linear semiflexible 
homopolymer chain confined in between two one dimensional stair shaped 
impenetrable surface (geometrical constraint) under good solvent condition 
in two dimensions. For two dimensional space, the surface is a line and 
polymer chain is constrained by the such surface. 

To analyze effect of the geometrical constraint on the conformational 
properties of the semiflexible polymer, we have chosen fully directed self 
avoiding walk model introduced by Privmann et. al [24, 25] and used gen- 
erating function technique to solve the model analytically. The results so 
obtained is used to discuss behaviour of the polymer chain in constrained 
geometry and also to compare the results obtained for conformational prop- 
erties of the polymer chain, when chain is in the bulk [26] and there is no 
constraint near it. 

The outline of the paper is as follows: In Sec. 2, we describe the lattice 
model of fully directed self avoiding walk and used it to model a semiflexible 
homopolymer chain confined in between constraints (stair shaped surface). 
We have solved the model analytically to calculate exact critical value of step 
fugacity required for polymerization of an infinitely long linear semiflexible 
polymer chain when the chain is confined in between the constraints. We 
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have discussed the variation of minimum critical value of step fugacity (for 
flexible chain) with separation between the constraints. We have also studied 
the variation of minimum and maximum possible values of the persistent 
length of the flexible and stiff chains respectively when it is confined in 
between the constraints. Finally, in Sec. 3, we discuss the results obtained. 



2 Model and method 

A model of fully directed self-avoiding walk [24, 25] on a square lattice 
has been used to calculate conformational properties a linear semiflexible 
homopolymer chain confined in between two impenetrable stair shaped sur- 
face under good solvent condition (as shown schematically in figure 1). The 
directed walk model is restrictive in the sense that the angle of bending has 
unique value, that is 90° (for present case) and directedness of the walk 
amounts to certain degree of stiffness in the walks of the chain because all 
directions of the space are not treated equally. However, directed self avoid- 
ing walk model can be solved analytically and therefore it gives exact value 
of the conformational properties of the polymer chain. Since, we consider 
fully directed self avoiding walk (FDSAW) model and therefore, walker is 
allowed to take steps along +x, and +y directions on a square lattice in 
between the constraints. 




Figure 1: This figure shows a walk of an infinitely long linear semiflexi- 
ble polymer chain confined in between two constraints (impenetrable stair 
shaped surface). All walks of the chain starts from a point O on the con- 
straint. We have shown three different cases viz. (A), (B) and (C) having 
separation (n) between the constraints along an axis three, four and five 
monomers (steps) respectively. The separation between the constraints have 
been defined on the basis of the fact that how many maximum number of 
steps a walker can successively move along any of the +x or +y direction. 
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The walks of the chain starts from a point O, located on an impenetrable 
surface and walker moves through out the space in between the two surface 
(as we have shown schematically in figure (1) that a walk of the polymer 
chain confined in between two surface for three different values of separation 
(n) between the two surface or the constraints). 

The stiffness of the chain is accounted by associating a Boltzmann weight 
with bending energy for each turn in the walk of the polymer chain. The 
stiffness weight is fc(= exp(—(3eb); where (3 = is inverse of the temper- 
ature, eb(> 0) is the energy associated with each bend in the walk of the 
chain, kb is Boltzmann constant and T is temperature) . For k = 1 or e& = 
the chain is said to be flexible and for < k < 1 or < e& < oo the polymer 
chain is said to be semiflexible. However, when 65 — > 00 or k — > 0, the chain 
has shape like a rigid rod. 

The partition function of a semiflexible polymer chain can be written as, 

^)=E^:r e (i) 

all walks of N steps 

where, Nb is the total number of bends in a walk of N steps (monomers), g 
is the step fugacity of each monomer of the chain. 

The partition function of an infinitely long linear semiflexible homopoly- 
mer chain confined in between the constraints (as shown schematically in 
figure 1A) can be calculated using the method of generating function tech- 
nique. The components (as shown in figure 2) of the partition function, 
Z^{k,g) (we have used here suffix three because in this case (i. e. figure 
1A) maximum step that a walker can move successively in one particular 
direction is three) of the chain can be written as, 



and 



Xi = g + kgY 3 
X 2 = g + g(X 1 + kY 2 ) 

X 3 =g + g(X 2 + kY 1 ) 

Y 1= g + kgX 3 
Y 2 =g + g(kX 2 +Y 1 ) 

Y 3 =g + g(kX 1+ Y 2 ) 



(2) 

(3) 
(4) 
(5) 
(6) 

(7) 



On solving Eqs. (2-7), we get the expression for X\(k,g) and Y 2 (k,g). 
In obtaining the expression for Xi(k,g) and Y 2 (k,g), we have solved a ma- 
trix of 2nX2n (n = 3, for present case i. e. figure 1A). Thus, we have 
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Figure 2: The components of the partition function is shown graphically in 
this figure. Term X m indicates sum of Boltzmann weight of all the walks 
having first step along +x direction, suffix m(l < m < n) indicates maxi- 
mum number of steps that walker can successively take along +x direction 
and n is separation between the two constraints. Similarly, we have defined 
Y m , where first step of the walker is along +y direction. In this figure, (B) 
and (C) graphically represents recursion relation for Eqs.(2&6) respectively. 



exact expression of the partition function for the semiflexible polymer chain 
confined between the constraints (as shown in figure 1A), can be written as, 

Z 3 (k,g) = X 1 {k,g) + Y 2 (k,g) = ~ , _ kg _ k 2 g2 _ kg3 + k s g s (8) 

From singularity of the partition function, Z^(k,g), we obtain critical value 
of step fugacity required for polymerization of an infinitely long linear semi- 
flexible polymer chain in between the constraints (as shown in figure 1(A)). 

The persistent length of the chain can be defined using tangent-tangent 
correlation function < t(s).t(s') >~ exp(—\s — s'\/l p ) [27]. The tangent 
vector t(s) can be defined as dr(s)/dt, here r(s) is parameterized in terms 
of arc length s of the polymer chain [27]. However, we have defined the 
persistent length l p , as the average length of the polymer chain between its 
two successive bends [26, 28]. Thus, the persistent length is calculated by 
us using following relation, 

C N > a _ dLogZsQ^g) dLogZ^g) 

h ~ <N b > ~ {9 dg )/[k dk j m 

(9) 

Where, a is lattice parameter of the square lattice and we have chosen its 
value unity for mathematical sake. 

The method discussed above can be used for different values of spacing 
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(n) between the constraints and size of the matrix needed to solve in calculat- 
ing partition of the chain confined in between the constraints is 2nX2n. We 
have calculated exact expression of the partition function for n (3 < n < 19). 

We have plotted the maximum value of the persistent length (l p (Max.) = 
l p (k = kMin.,9c(k = kMin.)) accessible to the stiff polymer chain confined 
in between the constraints for different values of n (3 < n < 19) in figure 
(3) and found that maximum value of the persistent length scales linearly 
with n as, l p = 0.66272 + 0.40723n. The value of kMin. is determined from 
singularity of the partition function Zz(k,g) and it gives least possible value 
of stiffness weight or maximum value of bending energy of the chain for 
which an infinitely long linear semiflexible chain is polymerized in between 
the constraints for a given value of n. Since, stiffness weight of the chain 
is related with bending energy (e&) as, k = exp(—f3et), therefore, we have 
calculated the maximum value of bending energy accessible to the stiff chain 
confined in between the constraints from minimum value of k{= kMin) and 
plotted maximum value of bendine enererv with n in figrure f3l. 




n n 

Figure 3: The exact value of the persistent length (l p ) is shown for different 
values of n (the maximum number of steps that a walker can successively 
move along either +x or +y direction in between the constraints). We have 
found that maximum value of the persistent length scales linearly with n(3 < 
n < 19). The maximum value of the persistent length (l p (Max.) = l p (k = 
kMin.,9c(k = kMin.)) follows a scaling relation, l p = 0.66272 + 0.40723n. 
We have also calculated the maximum value of the bending energy (/3e^ = 
—ln(k = kMin.)) associated with an infinitely long linear semiflexible chain 
which can be polymerized in between two constraints for different values 
of n. The maximum value of the persistent length and bending energy 
corresponds to a possible stiffer chain can be polymerized in between the 
constraints for given values of n 

However, minimum value of the persistent length accessible to the con- 
fined chain and minimum critical value of the step fugacity accessible to the 
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chain is shown for various values of n (3 < n < 19) in figure (4). 

2 r 
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Figure 4: In this figure, we have shown minimum value of the persistent 
length (Ip(Min) = l p (k = l,g c (k = 1)) and minimum critical value of step 
fugacity (g c (Min) = g c (k = 1)) required for polymerization of an infinitely 
long linear semiflexible polymer chain in between two constraints for differ- 
ent values of spacings (n) between the constraints. In other words, we have 
shown critical value of step fugacity and persistent length of a flexible poly- 
mer chain can be polymerized in between the constraints for n(3 < n < 19). 



3 Result and discussion 

We have considered an infinitely long linear semiflexible homopolymer chain 
confined in between two impenetrable stair shaped surface (constraint) in 
two dimensional space under good solvent condition. We have used fully 
directed self avoiding walk model to study effect of geometrical constraint 
imposed on the polymer chain and solved the model analytically to calculate 
exact expression of the partition function for few values of separation (n) 
between the constraints 3 < n < 19. We have found that maximum value 
of the persistent length scales linearly (l p = c + m * n, where c=0. 67745, 
m=0. 40472 and n — > oo). Thus, the persistent length of a possible stiffer 
chain which can be polymerized in between the constraints scales linearly 
with spacing (n) between the constraints. We have estimated value of c and 
m for spacing (n — > oo) between constraints using linear extrapolation (as 
shown in figure 5). We expect that exact maximum value of the persistent 
length will increase linearly for all values of n(> 3). 

We have also analyzed variation of the maximum value of the bending 
energy of the polymer chain confined in between the constraints and found 
that it increases with the increase of separation between the constraints 
for stiffer chains. The minimum value of the persistent length of a flexible 
polymer chain varies with separation between the constraints and appears 
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Figure 5: We have shown in this figure different value of c and m for few 
value of n and used linear extrapolation to obtain value of c and m for 

n — > oo. 

to approach a value less than 2 for n — >■ oo. However, critical value of step 
fugacity (g c (Min.)=g c (k = 1)) required for polymerization of an infinitely 
long linear flexible polymer chain in between the constraints approaches to 
0.5 as n — > oo. The minimum critical value of step fugacity is 0.5 (i. e. 
for k = 1) for a long linear polymer chain when chain is modeled as a fully 
directed self avoiding walk on a two dimensional square lattice in absence of 
the constraints, i. e. g c (k) = jt^ [26]. Since, g c {Min.) is critical value of 
step fugacity which is required for polymerization of an infinitely long linear 
flexible polymer chain on a two dimensional square lattice for fully directed 
self avoiding walk model and whence remain un-affected due to geometrical 
constrain imposed by stair shaped surface as n — >■ oo. 

However, minimum value of the persistent length is less than 2 for n — > oo 
while the persistent length (l p = 1 + k~ l ) [26] has value 2 for a flexible 
polymer chain (i. e. for k = 1) when there is no constraints near the chain. 
The persistent length is a local property, therefore, its value depends on the 
type of lattice (square, rectangular etc.) chosen to model the chain, space 
dimensionality and on the fact that polymer chain is ideal or self avoiding 
type. We have found that the value of l p (Min.) is less than 2 for n —> oo 
and this value is different from 2, when there is no constraints near the chain 
because for n — > oo there is still one constraint near the chain and confining 
to it or not allowing it to occupy space below the constraint at which the 
chain is grafted. 
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